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Abstract 

\ From the deep-inelastic momentum sum rule and the trace anomaly of the 

' energy-momentum tensor, I derive a separation of the nucleon mass into the 

contributions of the quark and gluon kinetic and potential energies, the quark 
. masses, and the trace anomaly. 
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The nucleoli derives its mass (939 MeV) from the quark-gluon dynamics of its underlying 
structure. However, due to the complexity of the low-energy Quantum Chromodynamics 
(QCD), a more detailed understanding of the nucleon mass seems difficult. The lattice QCD 
is successful in reproducing the measured mass from the fundamental lagrangian but the 
approach provides little insight on how the number is partitioned between the nucleon's 
quark and gluon content. Years after the advent of QCD, our knowledge of the nucleon's 
mass structure mostly comes from models: non-relativistic quark models. Bag models, the 
Skyrme model, string models, the Nambu-Jona-Lasinio model, to just name a few. Though 
all the models are made to fit the mass of the nucleon, they differ considerably on the account 
of its origin. Depending on different facets of QCD the models are created to emphasize, 
the interpretations of the nucleon mass often go opposite extremes. 

In this Letter I show that an insight on the mass structure of the nucleon can be pro- 
duced within QCD with the help of the deep-inelastic momentum sum rule and the trace 
anomaly. The result is a separation of the nucleon mass into the contributions from the 
quark, antiquark, gluon kinetic and potential energies, the quark masses, the gluon trace 
anomaly. Numerically, the only large uncertainty is the size of {P\msSs\P) , the strange 
scalar charge of the nucleon. Some implications of this break-up of the masses are discussed 
following the result. 

Let me begin with the energy- momentum tensor of QCD, 

T^"" = ^^W^^Y'^tp + ^g^T^ - F^T"^, (1) 

where ijj is the quark field with color, fiavor, and Dirac indices; F^'^ is the gluon field strength 
with color indices and F"^ = and all implicit indices are summed over. The covariant 

derivative D = D — D , with D = d + igA^^ and D = d — igAi^, where A^^ = A^f^ is 
the gluon potential. The symmetrization of the indices and v in the first term is indicated 
by (/iz/). Eq. (|I]) is quite formal, for it contains neither the gauge fixing and ghost terms, nor 
the trace anomaly. The first type of terms are BRST-exact and have vanishing physical 
matrix elements according to the Joglekar-Lee theorems I will add the trace anomaly 
explicitly when the renormalization issue is dealt with. 

A few results about the energy-momentum tensor are well-known. First of all, it is a 
symmetric and conserved tensor, 

rpt^y = T'^'^, d^T^"" = 0. (2) 

Because of the second property, the tensor is a finite operator and does not need an overall 
renormalization. All the fields and couplings in Eq. (|l]) are bare and their divergence are 
cancelled by the standard set of renormalization constants. The only complication is the 
gluon part of the tensor cannot be renormalized with a vanishing trace (see below). Second, 
the tensor defines the hamiltonian operator of QCD, 

^^QCD = / f/'xr°0(0,f), (3) 

which is also finite and scale- independent. Third, the matrix element of the tensor operator 
in the nucleon state is 

(^P\T^"'\P) = pt^P^/M, (4) 
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where \P) is the nucleoli state with momentum P^^ and is normahzed according to {P\P) = 
{E /M){2n)^S^{0), and E and M is the energy and mass of the nucleon, respectively. Lastly, 
the trace of the tensor is 

^ 1^,. [(^ ^ ^^^^^^ ^ ^ (5) 

where m is a quark mass matrix, 7^ is the anomalous dimension of the mass operator, 
and P{g) is the /3-function of QCD. At the leading order P{g) = — /?o5'^/(47r)^ and Pq = 
(11 — 2nj/3), where Uf is the number of flavors. The second term is called the trace anomaly 
and is generated in the process of renormalization. 
According to the above, the mass of the nucleon is 

_ {P\Jd^xT^%0,x)\P) , 

M = ^ ^^^^^^ ^' ^ = (T"") , (6) 

in the nucleon's rest frame. Although I formally work with the matrix elements of the 
nucleon, it actually is the difference of the nucleon matrix elements and the vacuum matrix 
elements that enters all the formula (the vacuum has zero measurable energy density). 
According to (^), a mass separation can be found through a decomposition of T'^'^ into 
various parts, which are then evaluated with the deep-inelastic momentum sum rule and the 
scalar charge of the nucleon. [Note that the parts of the energy-momentum tensor are not 
separately conserved, so the breaking of the nucleon energy cannot be Lorentz covariant.] 
First of all, let me decompose the T^'^ into traceless and trace parts, 

where T^'^ is traceless. According to Eq. (|), I have, 

(^P\f^^^\P) = (^pt^pi^ _ ^M^g^''')/M, (8) 

(PlT'^n^) = \9^''M. (9) 
Combining Eq. (P) with the above three equations, I get, 

(f = \m, (10) 

(f = \m. (11) 

Thus 3/4 of the nucleon mass comes from the traceless part of the energy-momentum tensor 
and 1/4 from the trace part. The magic number 4 is just the space-time dimension. This 
decomposition, a bit like the virial theorem, is valid for any bound states in field theory! 

The traceless part of the energy-momentum tensor can be decomposed into the contri- 
bution from the quark and gluon parts, 

ff^u^f^u^f^u^ (12) 

where 
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fl^'' = ^^zD^'y^^ - Ig'^^m^, (13) 

f^u ^ ^gt^^p'i _ (14) 
y 4 

Although the sum of Tj^'^ and T^'^ with bare fields and bare couplings is finite (now neglecting 
the trace anomaly), the individual operators are divergent and must be renormalized. Under 
renormalization, they mix with each other and with other BRST-exact and the equations 
of motion operators, which have vanishing physical matrix elements 0J^. For my purpose, 
I regard both operators renormalized and dependent on a renormalization scale /i^. Define 
their matrix elements in the nucleon state, 

{P\fg^^\P) = a{^?){P^'P'' - ^g^'^M^)/M, (15) 

l^pm^\P) = (1 _ a{^?)){P>'P'' - -g^"'M^)/M, (16) 
y 4 

where I have used Eq. (|^) to get the second equation. The constant a(/i^) is related to the 
deep inelastic sum rule 0, 

«(/^^) = / 2;[g/(a;, /x^) + g/(x, /i^)]rfx, (17) 

where the sum is over all quark fiavors and g/(x,/i^) and g/(x,/i^) are quark momentum 
distributions inside the nucleon in the infinite momentum frame or light-front coordinate 
system. Again, according the Eq. (^), I find the contribution to the nucleon mass, 

(T^) = (18) 

(f°°) = ^(l-a(/.2))M. (19) 

Finally, I turn to the trace part of the energy-momentum tensor T^'^ . According to 
Eq. (|^), I decompose it into Tjf^ and T^'^, the mass term and trace anomaly term, respec- 
tively. Both operators are finite and scale independent. If I define, 

b = 4(f oo)/M, (20) 

then according to Eq. (0), the anomaly part contributes, 

{f^') = \{l-b)M. (21) 

Thus, the energy-momentum tensor T^'^ can be separated into four gauge-invariant parts, 
T^", f^", f^^ and f^". They contribute, respectively, 3a/ A, 3(1 - a)/4, 6/4, and (1 - 6)/4 
fractions of the nucleon mass. The corresponding breakdown for the hamiltonian is, Hqcb = 
H'^ + Hg + H'^ + Ha, with 
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^ipmip^ , (22) 

(23) 
(24) 
(25) 

where I have consistently neglected 7^ and the higher-order terms in l3{g). One can put 
them back if a higher precision analysis becomes necessary. I also have taken nj = 3. [Note 
that the heavy quarks do contribute to the mass term, the kinetic and potential energy term, 
and the trace anomaly term. However, the contributions cancel each other in the limit of 
rrif 00, and thus for simplicity I neglect them.] If I rearrange the mass terms by defining, 

Hg = J d^x ^pi-iB ■ a)^p, (26) 
H^= [ d^x i/j7m{j, (27) 



then the QCD hamiltonian becomes, 

-f^QCD = Hq + Hm + Hg + Ha- (28) 

Here Hq (Eq. (^61) ) represents the quark and antiquark kinetic and potential energies and 
contributes 3(a — 6)/4 fraction of the nucleon mass. Hm (Eq. (pT])) is the quark mass term 
and contributes b fraction of the mass. Hg (Eq. (^)) is the normal part of the gluon energy 
and contributes 3(1 — a)/4 fraction of the mass. Finally, Ha (Eq. (^)) is the gluon energy 
from the trace anomaly. It contributes (1 — 6)/4 fraction of the mass. 

To determine the decomposition numerically, I need the matrix elements a and b. The 
deep- inelastic scattering experiments have determined a(/i^) with an accuracy of a few per- 
cent. Using a recent fit to the quark distributions 0, 

aj;js(lGeV^) = 0.55 (29) 

where MS refers to the modified minimal subtraction scheme. 
Without the heavy quarks, the matrix element b is, 

bM = {P\muuu + mddd\P) + (P|m,ss|P) (30) 

The first term is the vrA^ a-term apart from a small isospin-violating contribution of order 2 
MeV. A most recent analysis gave a magnitude of 45±5 MeV for this term P]. So the only 
unknown in our analysis is the strange scalar charge {P\msSs\P) in the nucleon. There are 
model calculations for this quantity in the literature 0. Here I choose to estimate it using 
two standard approaches, though both of them are not completely satisfactory. 

In the first approach, the strange quark mass is considered small in the QCD scale and 
so the chiral perturbation theory can be used to calculate the SU(3) symmetry breaking 
effects. A recent second-order analysis on the spectra of the baryon octet combined with 
the measured a-term yields 
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(P|SS|P) ~ 0.11 X (P|MM + rfrf|P) (31) 

~ 0.77, (32) 



where in the second hne, I have used (m„ + md)/2 ~ 7 MeV at the scale of 1 GeV^ [jlO 
Taking the strange quark mass to be 150 MeV at the same scale, I get. 



Using Eq. (pT]), I have. 



hM ~ 160 MeV, (33) 



(P|— p2|P) = -693 MeV. (34) 



In the second approach, the strange quark is considered heavy in the QCD scale. Using 
heavy-quark expansion, it was found [pJ]| , 

{P\mQQQ\P) = -^{P\'^F'\P). (35) 



Thus the strange quark contribution in 



{P\i)mi) + ^I^P^IP) = M, (36) 

which is an explicit form of Eq. (|ll|), cancels. From the above equation and the cr-term, I 
find, 

rv 

(P|— p2|p) = -740 MeV. (37) 

TT 

This yields a strange matrix element {P\msSs\P) = 62 MeV. Together with the cr-term, I 
determine, 

bM = 107 MeV. (38) 

The complete result of the mass decomposition at the scale of /i^ = IGeV^, together with 
the two numerical estimates, is shown Table 1. I have not shown the errors due to omission 
of higher-order perturbative effects and errors on the a-term and current quark masses. 
The total effect on individual numbers is about 5 to 10 MeV. Thus I have rounded up the 
numbers to nearest 10 MeV. The largest uncertainty is from the matrix element {P\msSs\P) , 
which could be larger than the difference of the two estimates shown. Nevertheless, I will 
argue below that the total strange contribution to the nucleon mass is quite small and with 
a smaller uncertainty. 

The following comments can be made with regard to the numerical result. 

• The quark kinetic and potential energies contribute about 1/3 of the nucleon mass. 
Because the quark kinetic energy must be very large when confined within a radius of 
1 fm, there must exists a large cancellation between the kinetic and potential energies. 
This may not be entirely surprising in the presence of strong interactions between 
quarks and gluons. Such strong interactions are clearly at the origin of the chiral 
symmetry breaking, modeled, for instance, by the Nambu-Jona-Lasinio [H. 
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• The decomposition of the quark energy into different flavors is possible. Taking the 
number 270 MeV (the rris limit) as an example, I find the up-quark energy 
in the proton is 250 MeV using the momentum fraction carried by up quark 0.375 
[0], the down quark energy, 105 MeV, and strange quark energy, —85 MeV. Further 
decomposition into valence and sea contribution cannot be made without knowledge 
of the separate valence and sea contributions to the scalar charge. 

• The quark mass term accounts for about 1/8 of the nucleon mass. About half of 
which or more is carried by the strange quark. The contributions from the up and 
down quarks are well determined by the a-term. 

• The normal gluon energy is about 1/3 of the nucleon mass and the trace anomaly part 
contributes about 1/4. From these two, I deduce the color-electric and color-magnetic 
fields in the nucleon separately (take as(lGeV)~ 0.4), 

(P|EV) = 1700 MeV, (39) 
(P|BV) = -1050 MeV. (40) 



So the magnetic-field energy is negative in the nucleon! This of course is due to a 
cancellation between the quark's magnetic field and that of the vacuum. The electric 
field in the vacuum is presumably small, however, it is large and positive in the nucleon. 
This behavior of the color fields in presence of quarks is very interesting, it may help 
to unravel the structure of the QCD vacuum. 

In the chiral limit, the gluon energy from the trace anomaly (M/4) corresponds exactly 
to the vacuum energy in the MIT bag model ||T^ . The role of such energy in the model 



is to confine quarks. Thus we see here a clear way through which the scale symmetry 
breaking leads to quark confinement. To keep the confinement mechanism, a model 
must include Ha, the anomaly part of the hamiltonian. 

• The strange quark contributes about —60 MeV through trace anomaly. When adding 
together with the kinetic and potential energy contribution —85 MeV and the mass 
term 115 MeV (the limit) the total strange contribution to the nucleon mass 

is a mere —30 MeV. (The other limit gives a total of —45 MeV.) The smallness of the 
contribution is, to a large extent, insensitive to the matrix element {P\msSs\P) . 

To summarize, we have find a separation of the nucleon mass into contributions from 
the quark kinetic and potential energy, gluon energy, and the trace anomaly. The largest 
uncertainty is from the strange matrix element {P\msSs\P). The result has interesting 
implications on the quark-gluon structure of the nucleon and on the response of the QCD 
vacuum to color charges. 

I wish to thank D. Freedman, F. Low, and K. Johnson for useful discussions and sugges- 
tions. 
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TABLES 



TABLE L A decomposition of the nucleon mass into different contributions. The matrix ele- 
ments a and b are defined in Eqs. (p!5|) and (pO[). 



mass type 
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quark energy 
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'ipmtp 


b 
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gluon energy 
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